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Abstract. We report the serendipitous discovery of an 
extremely red object in the Hubble Deep Field South 
(HDFS) Test NICMOS (Near Infrared Camera and Multi 
Object Spectrograph) field of view. The object is re¬ 
solved in the NICMOS image and has a light profile 
very well described by an r 1//4 law with effective radius 
re = O'/20 ± 0.05 and Hab = 21.7 ± 0.1 magnitudes. In 
contrast, the galaxy is undetected in the R and I band 
ground based images taken at the CTIO 4 m Blanco Tele¬ 
scope, giving a lower limit to the color of (R-H)ab> 3.9 
and (I-H)ab> 3.5 at the 95 % confidence level. The col¬ 
ors of a range of synthetic galactic spectra are computed, 
showing that the object is likely to be an “old” elliptical 
galaxy at redshift z <; 1.7. Alternatively the colors can 
be reproduced by an “old” elliptical galaxy at somewhat 
lower redshift (<; 1) with significant amount of dust, or by 
a younger galaxy at higher redshift. This object represents 
a very interesting target for future VLT observations. 

Key words: Galaxies: elliptical and lenticular, cD- 
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1. Introduction 

The joint effort of the Hubble Space Telescope (HST) 
and large ground based telescopes in the last few years 
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has produced a real breakthrough in our understanding 
of the history of formation and evolution of galaxies. 
Deep photometric multicolor imaging such as the Hub¬ 
ble Deep Field (HDF, Williams et al. 1996) together 
with spectroscopic information (e.g., Lilly et al. 1995; 
[Steidel et al. 199(| [Lowenthal et al. 1997 ) have started to 
sketch a sort of cosmic timetable that provides impor¬ 
tant constraints to theories for the formation of cosmic 
structure (Steidel et al. 1998; Baugh et al. 1998). Besides 
global properties ( Madau et al. 1996Q , the superb resolu¬ 
tion of the HST allows us to study the evolution of differ¬ 
ent morphological types separately, in particular to clarify 
the history of elliptical galaxies (E/SO). The intermediate 
redshift spectroscopic data (e.g. van Dokkum et al. 1998) 
suggest a mostly passive evolution of old stellar popula¬ 
tions that seems t o be confirmed by photometric studie s 
at higher redshift ( [Schade et al. 1997 ; Driver et al. 1998 ). 

To push the investigation to redshift significantly 
greater than 1, infrared (IR) photometry (corresponding 
to optical rest-frame emission) is needed (e.g. Maoz 1997; 
[Dickinson 1998| ). The extremely red objects with (R—K)> 

5 (equivalently (R— K)ab > 3.3) found with IR sur¬ 
veys are generally thought to be high redshift elliptical 

galaxies (e.g., Hu fe Ridgway 1994 ; Spinrad et al. 1997 ). 

However, due to their subarcsecond sizes, HST is needed 

to explore in detail the morphology of these objects 
(Graham & Dey 1996). 

In this Letter we report the discovery of a “R-H 

dropout” (Figs. |l] a, b and c in the HDFS test image. 
The morphology, the size and the colors (R-H and I-H) 
strongly suggest the object to be an “old” elliptical galaxy 
(or bulge-dominated object) at z 1.7 and therefore with 
a very high formation redshift (zj), somehow similar to the 
“old” red galaxy found by Spinrad et al. (1997) at redshift 

z = 1.55. Many extremely red objects have been found in 


the last years (see e.g. [Elston et al. 1988[ [Dey et al. 1995 
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Fig. 1. a. HST NICMOS F160W “drizzled” image (7680 
s). b. CTIO 4m R (3000 s). c. CTIO 4m I (2800 s). North 
is up and East is left in the CTIO exposures, while the 
NICMOS image is slightly rotated. The frames are « 35" x 
45". The boxes are centered on the position of the galaxy. 




Moustakas et al. 1997 ; Stanford et al. 1997 ) using ground 
based IR photometry, thus with poor morphological infor¬ 
mation. At the opposite, the object we present in this Let¬ 
ter, thanks to the HST angular resolution, is resolved and 
we can measure the light profile and the effective radius 
(similar objects are being detected with NICMOS, e.g. 
Me Carthy et al, 1998, but the luminosity profile has not 
been measured). The luminosity profile is well described 
by an r 1 / 4 law. Models with dust are also considered to ex¬ 
plain the extremely red colors, given the growing evidence 
for significant amounts of dust in high redshift galaxies 
( [Softer et al. 1998 ; Fimatti et al. 1998| ). The Hubble con¬ 

stant is assumed to be 100 h km/s/Mpc, with h = 0.65 
where needed. 


2. Photometry 

During the HDFS test program (Williams et al. 1997) 
images of the candidate HDFS held were taken with the 
F160W filter of the Near Infrared Camera and Multi Ob¬ 
ject Spectrometer (NICMOS) on board the HST (Fig. |lj 
a), and in the optical R band (Fig. [l] b) with the 4m 
Blanco Telescope at the Cerro Tololo Inter-American Ob¬ 
servatory (CTIO). Follow up I band photometry was ob¬ 
tained on May 13 1998 at the same telescope (Fig. [i] c). 

2.1. The Infrared NICMOS Imaging 

Twelve NICMOS Camera 3 (NIC3) images of the held 
were taken in MULTIACCUM mode (640 s each). The 
images were obtained at six different pointings in order to 
correct for bad pixels and to recover the information lost 
to undersampling of NIC3 (approximately 0"2 pixel size). 

The “pedestal”, an unpredictable bias that affects the 
NICMOS images (see the NICMOS Image Anomalies web 
page at URL www.stsci.edu/ftp/instrument_news/ NIC- 
MOS/nicmos_anomalies.html), was removed using the 
Pedestal Estimation Software by R. P. van der Marel, 
and scripts developed by one of us (MED). The dithered 
images were combined on a subsampled grid (O'.'l pixel 
size, shown in Fig. |l] a) using the DRIZZLE iraf/stsdas 


task (Fruchter & Hook 1998). A bad pixel and cosmic ray 
mask was obtained by using the iraf/stsdas package 
ditherII. 


The object appears resolved on the hnal drizzled im¬ 
ages, but to be sure that this was not an artifact of the 
reduction procedure we performed a y 2 * test on the sin¬ 
gle calibrated frames. For each image we measured the 
FWHM along the x and y axis of the galaxy and the two 
nearby stars (using the MIDAS command center/Gauss). 
We estimated the error on the FWHM to be the standard 
deviation of the values measured on the different frames 
(about 0"04). The reduced y 2 of the values with respect 
to their average is therefore 1 (by definition). The reduced 
y 2 of the measured FWHMs computed with respect to the 
average stellar FWHM is y 2 = 11.1, showing that the ob¬ 
ject is resolved. 

The luminosity profile of the galaxy (Fig. |j) was ob¬ 
tained by fitting elliptical isophotes to the image with a 
modified version of the MIDAS command fit/ell3 specifi¬ 
cally designed to deal with undersampled images (Mpller, 
Stiavelli & Zeilinger, 1995). 

In order to fit a luminosity profile law to the data point 
it is crucial to determine the PSF very accurately. An 
unexpected deformation of the NICMOS dewar pushed 
NIC3 out of the range reachable by the focusing device 
(Pupil Alignment Mechanism) and therefore NIC3 is al¬ 
ways somewhat out of focus. Thermal breathing of the 
instrument (that is significant during long exposures) and 
the combination of multiple dithered exposures add fur¬ 
ther uncertainty and broadening to the PSF. For these 
reasons the stars in the field are a better approximation 
of the real PSF than the PSFs obtained with Tiny Tim 4.4; 
the main residual difference may be caused by position- 
dependent features of the NIC3 PSF. In order to have an 
estimate of the errors due to the PSF, we fitted the profile 
with different PSFs: the two nearest stars and two “syn¬ 
thetic” ones computed as follows. A PSF was calculated 
using Tiny Tim 4.4 (using a 15 mas jitter) on a subsam¬ 
pled grid (0"02 pixel size) in the same position of a star 
in the HDFS test image. It was then rebinned to the scale 
of the drizzled image in order to obtain a PSF centered in 
the same sub-pixel of the real star. Under the assumption 
that the “true” PSF (S) is the convolution of the Tiny Tim 
PSF (T) with a position independent broadening function 
( B ) due mostly to the breathing of NICMOS and the driz¬ 
zling of the single dithered images, we derived the func¬ 
tion B using the MIDAS fourier transform commands. The 
broadening functions B\ and B 2 were derived for the two 
nearest stars. Using a Tiny Tim PSF (Tg) centered on the 
galaxy center and the broadening functions we obtained 
two “synthetic stars” synl=Xg o Bi and syn2= Tg o H 2 . 
As we were interested in checking the robustness of the re¬ 
sults with respect to PSF uncertainties, we performed the 
fits using both the two real stars and the two “synthetic” 
stars. 
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Fig. 2. IR (F160W) luminosity profile of the galaxy 
(crosses). The PSF-convolved best-fit r 1 / 4 (a) and expo¬ 
nential (b) profiles are shown. Only a few error bars are 
shown for clarity. Four different fits are obtained by using 
different PSF profiles as described in the text. 

We fitted an r 1 / 4 law to the isophotal profile, using a 
code described, e.g., by Carollo et al. (1997). The fit ob¬ 
tained with the four PSFs are in excellent mutual agree¬ 
ment and reproduce the data very well (Fig. 0 a). If an 
exponential profile component is added, the best fit is 
obtained for a very small exponential component and a 
practically unchanged r 1 / 4 component. At the opposite, 
an exponential disk by itself gives a poor description of 
the profile (Fig. 0 b). The photometric parameters were 
also determined with a two-dimensional fit code, described 
in Treu et al. (1998). The total error is mainly due to 
residual PSF mismatches, to the specific technique (pro¬ 
files or 2D fit), and to errors in the sky subtraction. The 
PSF and the technique dependent uncertainties can be 
estimated from the scatter of all the photometric param¬ 
eters derived, while the sky-subtraction error has been 
estimated by varying the sky level (« 15 % in re and 
« 0.1 mag in m). The average values of the effective ra¬ 
dius and magnitude, obtained with different techniques, 
models and PSFs, are respectively re = 0'.'20 ± 0.05 and 
m=21.7 ± 0.1. 

2.2. The R and I ground-based photometry 

Ten dithered R band exposures (300 s each) were obtained 
with the Prime Focus CCD (PFCCD) Imager at the CTIO 
Blanco 4m telescope on May 2 1997, using the detector 
SITe 2K #6. The night was photometric with seeing 1"2 
. The same field was observed in the I band on May 13 
1998 (unfortunately almost Full Moon was present) with 
the Tyson-Bernstein Mosaic Imager. CCD #3, also a SITe 
2K, was pointed at the PFCCD field, obtaining 14 x 200 s 
dithered exposures. The night was photometric, with l'.'O 
seeing. The images were reduced and combined in stan¬ 
dard way using iraf. 

The photometric calibration of the R band was done 
by bootstrapping the photometry from a 0.9m to the 4m 
image. Three relatively bright (I ~ 18) stars near the NIC- 


MOS position were calibrated on the I-band frame by ref¬ 
erence to seven stars in Landolt (1992) SA 107. Both were 
also checked with the nominal zeropoints from the instru¬ 
ment manuals. The accuracy of the calibration is about 
0.1 mag, which is good enough for our purposes. 

The galaxy, which is quite evident on the NICMOS 
images, is undetected in both the R and I CTIO image 
(Figs. 0b and c) . Thus we can only set an upper limit to 
the R and I band luminosities. The nearest detected source 
is at ~ 3" on the R band image. The limiting magnitudes, 
computed by measuring the sky noise (JK) and correcting 
the resulting flux for the seeing losses, are iur > 25.9 and 
mj >24.7 AB mag at the two-sigma level. 

A more accurate limit can be given by considering the 
statistical distribution of the counts on the detector for 
a given source of intensity Iq and a given background of 
zero average and variance JK 2 . The 95 % confidence level 
(CL) limit is mq > 25.6 and nq >25.2 ABmag. The limits 
on the color are therefore (R-H)ab> 4-3 and (I-H)ab>3, 
at the two-sigma level, considering only the sky variance, 
but (R-H)ab>3-9 and (I-H)ab>3-5, at the 95 % CL with 
the more accurate algorithm. 

3. Identification 

The morphology and the IR light profile suggest that the 
object is a high redshift “old” elliptical galaxy. But the 
strongest clues in favor of a high redshift elliptical galaxy 
are the colors. To better constrain this identification we 
have computed the (R-H)ab and (I-H)ab colors of a 
set of synthetic spectra of single burst elliptical galax¬ 
ies (Bruzual & Chariot, 1993; GISSEL 96 version) as a 
function of redshift. A standard Salpeter IMF and metal- 
licities Z = Z & and Z = 0.2 Zq have been used. The lower 
metallicity has been selected to be representative of the 
low metallicity environment that one could expect to find 
in primordial stellar populations. In Figs. 0 a and b the 
contour plots of (R-H)ab and (I-H)ab (isochromes) as a 
function of redshift and age of the synthetic galaxy are 
shown. For each given redshift the galaxy age has to be 
less than the age of the Universe, which is overplotted as 
thin lines for two different values of the cosmological pa¬ 
rameters fl and Da- In practice, only the points above 
the maximum of the lower branches of the two isochromes 
(I and R), and below the thin line for the selected cos¬ 
mology, are consistent with the observed colors. As can 
be seen, the colors support the identification as an “old” 
elliptical galaxy at redshift ^ 1.7 for the solar metallicity 
models. The lower limit in redshift is even more stringent 
if we consider low metallicity spectra, because the Balmer 
Jump is intrinsically shallower. A redshift z ^ 1.7 would 
also give an effective radius in the usual range for elliptical 
galaxies (^ 1 kpc). The colors are well reproduced also for 
much higher redshift (~7 15). 

A significant content of dust (which can be found in 
high redshift star-forming galaxies, see e.g. Soifer et al. 
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of galaxies have been photometrically identified in the 
vicinity of QSOs (Francis et al. 1996; Steidel et al. 1998). 
Therefore this galaxy may possibly be a companion to 
the QSO J2233-60, which is located « 8', i.e. « 2/i -1 
Mpc (physical distance) from the HDFS NIC3 field. This 
identification would then lead to 2 ~ 2.22 compatible 
with the colors and size that we have observed. At this 
redshift, Hab=21.7 would imply an absolute magnitude 
M v = -22.7 (fi=l, fl A = 0) and My = -23.3 (fi=0.3, 
Oy=0.7) which are normal for a 1.5 Gyr-old elliptical 1 !. 

Further spectroscopic and photometric investigation 
are needed to measure the redshift, to study in detail 
the stellar population, to clarify the role of dust, and 
to identify possible high redshift companions. The faint¬ 
ness of the object though, makes it extremely hard to 
measure the redshift with 4m class telescopes, while it 
will be feasible in a few hours with the ISAAC spec¬ 
trograph at the VLT (see the VLT-ISAAC web page 
www.eso.org/instruments/isaac). 


AGE (Gyrs) 

Fig. 3. a. Contour plot (thick lines) of the (R-H)ab 
(solid) and (I-H)ab (dashed) colors (isochrome) as a func¬ 
tion of redshift and age of an elliptical galaxy, computed 
using synthetic spectra of solar metallicity. The age of the 
Universe as a function of redshift is overplotted as thin 
lines for two different cosmological models. See text for 
discussion and details, b. As in a but the models have 
metallicity Z = 0.2 Zq. c. As in a but the spectra have 
been reddened for dust (Ay = 0.5, see text), d. As in c 
but the models have metallicity Z = 0.2 Zq. 


1998, Cimatti et al. 1998) could redden the colors of a 
given galactic spectrum and redshift. For this reason, the 
same isochromes were computed after reddening the spec¬ 
tra for dust absorption following the extinction law given 
by Cardelli et al. (1989) with Ay = 0.5, Ay = 3.1E(B-V). 
In Figs. H (c) and (d) the isochromes are shown in the 
low redshift (0-3) range. The colors are now reproduced 
by lower redshift models (* ^ 1), and therefore the ob¬ 
ject could be a high redshift “old” elliptical galaxy with 
significant dust content. 


4. Conclusions 

The object, that has been selected as an “R-H dropout” 
in the test image of the HDFS, is clearly resolved on 
the NICMOS image and shows an r 1 * * / 4 law profile with 
an effective radius of 0"20. The colors (R-H)ab and (I- 
H)ab ar e well reproduced by synthetic models of “old” 
elliptical galaxies at redshift z <; 1.7 or at higher red¬ 
shift (up to ~ 15) by younger ones (Sect. ||), implying 
therefore a very high formation redshift. Alternatively 
it could be a lower redshift (z )> 1) elliptical galaxy 
with a significant amount of dust. High redshift clusters 
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